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■ Abstract 

<N : 

I ' We give a self-contained presentation and comparison of two different al- 

^ ■ gorithms to explicitly solve quantum many body models of indistinguishable 

On I particles moving on a circle and interacting with two-body potentials of 1/ sin^- 

I type. The first algorithm is due to Sutherland and well-known; the second one 

Tij- ■ is a limiting case of a novel algorithm to solve the elliptic generalization of the 

, Sutherland model. These two algorithms are different in several details. We 

Q I show that they are equivalent, i.e., they yield the same solution and are equally 

^ ' simple. 
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> ' 1 Introduction 

•i-H . 

^ • 

d ■ We recently presented a novel algorithm to solve the quantum version of the ellip- 
tic Calogero-Moser-Sutherland system |L^]. In the trigonometric limit, such an 



algorithm was discovered already about thirty years ago by Sutherland ||Sul| , [Su2 
Somewhat surprisingly, the former algorithm in that limit reduces to one which is dif- 
ferent from Sutherland's, even though it yields the same solution and is equally simple. 
The purpose of this paper is to give a detailed and self-contained comparison of these 
two algorithms, including a proof of their equivalence. 

The Sutherland model is defined by the differential operator 

^ = -E^ + 2A(A-1) 5^ Vix,-x,) (1) 

j=l i l<j<k<N 
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with —IT < Xj < IT, N = 2,3, X > 0, and 

Vir) = . } (2) 
4 sm (r/2) 

This differential operator defines a self-adjoint operator on the Hilbert space of square 
integrable functions on [— vr, vr]^, providing a quantum mechanical model for in- 
distinguishable particles moving on a circle of lengthQ 2tt and interacting with a two 
body potential proportional to V{r) where A determines the coupling strength. (To 
be precise: This model corresponds to a particularly nice self-adjoint extension of this 
differential operator which, for A > 1, corresponds to the Friedrich's extension [ RS|| .) 



To solve this model amounts to constructing a complete set of eigenfunctions and 
corresponding eigenvalues of this Hamiltonian. 

The starting point for Sutherland's algorithm is the following 
Fact 1 [ Sul |: The wave functiorf^ 



*o(x)= n ^i^k-x,)^ (3) 

i<j<k<N 

with 

i){r) = sin(r/2) (4) 
is the ground state of the Sutherland Hamiltonian, H'^q = -Eq^o- 

Exploiting this fact, Sutherland constructed all other eigenfunctions / using the 
following ansatz, 

/(x) = vl/o(x)<|.(x) (5) 

where $ are symmetric polynomials (i.e. non-negative powers) in the variables zj = 
exp(ixj) ||Su2|| . The symmetric polynomials thus obtained are the so-called Jack poly- 



nomials which have been studied extensively in the mathematics literature, see e.g. 



Our algorithm is based on the following 

^To ease notation, we set the length of space to 2tt from the start. Of course, an arbitrary length 
L > can be easily introduced by rescaling xj {2'K/L)xj, H {2'k/L)'^H, etc. 

^To fix the phase of unambiguously one can interpret sin(r/2)'^ as limejo sin(r/2 + ie)'*', for 
example. Anyway, the phase ambiguities associated with the exponentiated sines are irrelevant here. 
In Appendix B.l we will have to be more careful about similar phase ambiguities in the functions 
i^(x; y) defined below. 
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Fact 2 ||L1|| : The function 



^, . ]\i<j<k<N'4'{xk-XjY]\^^^^^^^'4){yj-ykY 

F(x;y) = — — = — — = , (6) 

ip{r) as in Eq. obeys the following identity, 

j=l i l<j<k<N 

with V{r-) as in Eq. 

Note that we can write this latter identity as 

i/(x)F(x;y) = if(y)F(x;y) (8) 

where H is the differential operator in Eq. (|I]) but acting on different arguments x 
and y, as indicated. The idea of our algorithm is to take the Fourier transform of Eq. 
(|^) with respect to the variables y, and this yields an identity allowing to construct 
eigenfunctions and the corresponding eigenvalues (Proposition 0). 

It is interesting to note that Fact 2 holds true in the elliptic well (in this 

case, ifjir) is a Jacobi Theta function 'diir/2) with nome q = exp(— /?/2) and V{r) is 
Weierstrass' elliptic function p{r) with periods 2n and i/3) ||L1|| , in contrast to Fact 1 



5u3|| . For the convenience of the reader, an elementary proof of Fact 2 (in the trigono- 
metric case) is given in Appendix A. (This proof uses Fact 1; a self-contained proof 
valid also in the elliptic case will be given in p^ .) 



The plan of the rest of this paper is as follows. In Section 2 we review the Suther- 
land algorithm | {Su2| , mainly to introduce our notation. Section 3 contains a detailed 



description of our algorithm. In the final Section 4 we give the arguments which prove 
that both algorithms are equivalent, despite of various differences. Lengthy proofs are 
deferred to two Appendices. 



2 Sutherland's algorithm 



We use H-^o = Eq^o with Bui 



Eo = -X'N{N^ - 1) 



(9) 



and make the ansatz / = $\l/o. With that the eigenvalue equation Hf = Ef becomes 



H'^ = E'^ with E' = E-Eq and 



(10) 



i<j<k<N 

One now determines the action of H' on symmetric polynomials 

AT 

P j=l 

where0 

ni > n2 > . . . > > 
and the sum is over all permutations P of {1, 2, ... , A^}. Using the identity | |Su2| | 



(12) 



(e^^ + e'y) ( — — ) = e''^^ + e'''^ + 2 V e' 

u=l 

for k > 0, one obtains 

nj— nj.— 1 

H'Sn = E'^Sn + A ^ (rij - Uk) ^ S^-u^^^ 

l<j<k<N u=l 

where we introduced the notation 

{Ejk)e = Sji - Ski, j,k,i = l,2,...,N 

and defined 



fc-i 



[{k-u)x+i^y] 



(13) 



TV 



N 



3=1 l<j<k<N 3=1 



n. 



(14) 



(15) 



We now introduce the notation 



(16) 



Note that Uj in Ref. |Su2|] corresponds to rijv+i-j here. 
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for non- negative integers fi^^, and observe that there is a natural order on the set of 
all these fi (which we can identify with Nq^'^ ^^^'^), 

H < if /Ujfc < /Uj-fc for all j < k. (17) 

It is obvious that H'Sn is a finite linear combination of symmetrized plane waves 5*11-^ 
with /J. > 0. We thus can make the following ansatz for the eigenf unctions of H', 

<^n = J2c^S^-^. (18) 

with 

= if (n — < (n — /x)^ for at least one j < k. (19) 

The latter condition shows that there are only a finite number of non-zero c^, i.e., the 
$n are polynomials. Then H'^^ = E'^^ implies E' = E'^ and the following recursion 
relations for the coefficients c^, 

nj-rik 

[EL - K-Jc^ = ^ E E - ^)^- - - t)^ + 2^]^M-E.- (20) 

(we used the fact that the functions are linearly independent). We can set Cq = 1 
(this fixes the normalization of the eigenfunctions) and then determine the other 
recursively, which is possible provided that there is no resonance, i.e., if E'^ — E'^_^^ is 
non-zero. This is the case: we shall prove at the end of this Section that 

E'n - = E -"^-fct^^ - B)^ - - t)!^ + ^""^ - + 2^(^ - j)] (21) 

which is manifestly positive and shows that resonances indeed do not occur. This 
completes the construction of eigenfunctions and eigenvalues of the Sutherland model: 
Note that the symmetrized plane waves provide a complete orthonormal basis of the 
corresponding non-interacting Hamiltonian (obtained by setting A = 0), and we have 
constructed eigenfunctions /„ = $n^o and corresponding eigenvalues En = E'^^ + Eq 
with Eq in Eq. (^, which are one-to-one to this free solution which is known to provide 
a complete basis. 

Remark 1 We can write (cf. Eq. (|T|)) E'^ = + |^(^ + 1 - 2j)]^ - with 

^ 1 

E'o = J2-X\N + l-2jf. (22) 
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It is easy to show that E'^ = Eq (cf. Eq. {^)), which is somewhat remarkable and imphes 
the following simple form of the eigenvalues, 

N 2 

^n = $^(ri, + A[i(iV+l)-j]) . (23) 
i=i 

The novel algorithm in the next Section will yield this simple form of the eigenvalues 
directly. 



For the convenience of the reader, we conclude this Section with a 
Proof of Eq. (^: 
Eq. ( [T3| ) implies 

K - K-^ = Yl^'^^-^'^ + 2^^' + A(iv + 1 - 2j)] 

with 

k>j k<j 

Using 

AtjOj = Y ~ ^k) (25) 

we get 

K - K-f. = Yl - z^^) + ^K- - ^k) + 2A(A; - j)] 
j<k 

which proves Eq. (|2l|). □ 



3 The novel algorithm 

This algorithm is based on the following Proposition which, roughly speaking, is ob- 
tained by taking the Fourier transform of the remarkable identity in Eq. (§) with respect 
to y. 

Proposition 1 Let H be as in Eqs. (^-(^. Then 

oo 

i7F(x; n) = E„F(x; n) - 7 ^ ^ z/F(x; n + z/E,,) (26) 

^<j<k<N u=l 
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where 

F(x;n) =P„(x)*o(x), n G 
with \l/o as in Eqs. and 



(27) 



Pn(x) = lim 



i-. (2vr) 



ni 



<j<k<N 



(2vr) 



X _ Q^iyj-yk)-ik-j)e \ 
^ (28) 



Ej-fc as in Eq. i^n m Eg. and 

7 = 2A(A-1). 



(29) 



(Proof in Appendix B.l.) 

Remark 2 To see that these functions are well-defined, we note that they can be 
written as 



K(x) 



d6 ^ni 



with integration paths Cj : = e'^-'e'^^, — tt < Uj < tt, where e > is arbitrary. 



We now show that this proposition provides a solution algorithm: Eq. (|26| ) implies 
that the action of H on the functions -F(x; n) is triangular, i.e., HF{x; n) is a linear 
combination of functions F(x; n + ^) with /x > 0. We thus can make the following 
ansatz for eigenfunctions. 



/n(x) = ^ a/^F(x; n + /x) 

M>0 



(30) 



and then Hf^ = Ef^ implies 
^F(x;n + ^) 

M>0 



En+fi — E 



^<j<k<N u=l 



uE 



<jk 



0. 



We thus see that we get a solution of -ff/n = -E/n is we set E = E^ and determine the 
coefficients by the following recursion relations. 



l<j<k<N v=l 



<jk 



(31) 



which has triangular structure: we can set Oq = 1 (this fixes the normahzation), and 
then the other can be determined recursively in terms of the a^/ where /x' < n, at 
least if there is no resonance, i.e., if -En+zx — does not vanish. This is true due to 
the following 

Lemma 1 

N 

^n+M-^n = 5^^?+ Yl '^f^Jklirij - Uk) + X{k - j)] (32) 
j=l l<j<k<N 

with fij in Eq. which is manifestly positive provided that Eq. ( ^T^ ) holds true. 

(Proof in Appendix B.2.) 

Moreover, the following Lemma shows that the /„ are in fact symmetric polynomi- 
als, i.e., a finite linear combination of the functions S-a in Eq. (|11]). 



Lemma 2 The functions Vn in Eq. ( ^Bj ) all are symmetric polynomials in the variables 
Zj = exp {ixj) which are non-zero only if 

nj + nj+i + ... + nN>0 Vj = 1, 2, . . . . (33) 

They can be written as 

Vn{^) = J2Pn,mSm{^) (34) 
m 

with S'ni(x) as in Eq. (jTTlj, and the coefficients are 

where the sum Y^" is over all non-negative integers /i^fc, uji^ restricted by the following 
2N equations, 

N j-l N N 

t=i 1=1 i=j+i i=i 

and mi > m2 > ■ ■ ■ > ttln > 0, implying in particular that there are only terms such 
that 

N N 

= Y^^- ^^^^ 

j=l 3=1 
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(Proof in Appendix B.3.) 

Indeed, this Lemma implies the sum in Eq. ( |30D has only a finite number of non- 
zero terms (since there is only a finite number of /J. such that n' = n + /x obeys all the 
conditions in Eq. (P^D), and thus the /„ are a finite number of terms each of which is 
a finite linear combination of functions S'n in Eq. (|ll]). 



4 Conclusions 



We can summarize our discussion in the previous Sections as follows. 

Proposition 2 For each n G such that ni > n2 > ■ ■ ■ > > 0, the standard 
algorithm reviewed in Sections 3 and the novel one presented in Section 4 both yield 
an eigenfunction /„ of the Sutherland Hamiltonian H in Eq. In both cases, this 
eigenfunction is of the form 

/„(x) = <l'„(x)vl/o(x) 

with in Eqs. and a symmetric polynomial in the variables Zj = exp (ixj), 

and the corresponding eigenvalues E^i are given in Eq. ^Bj). 



It thus follows from Theorem 3.1 in Ref. |^ that, for non-degenerate eigenvalues 
E^, the eigenf unctions /„ obtained with the two algorithms are equal (up to normal- 
ization), and the functions are proportional to the so-called Jack polynomials (see 
Section 2 in Ref. @ for details^). We feel that this is quite remarkable since, even 
though the two algorithms look somewhat similar and both yield the same solution, 
there are several differences in details: 



• The building blocks of the eigenf unctions in the novel algorithm are the functions 
Vn defined in Eq. (|28D and not the plane waves in Eq. (|TT1) . 

• With the standard algorithm, one obviously obtains eigenfunctions with polyno- 
mials $n which have the form 

m<n 

"'The latter Reference actually seems to suggest that this is true even for non-degenerate eigenvalues. 
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where the partial order here is defined as 



k k 

m < n :^ < ^rij Vfc = 1, 2, . . . (38) 

j=i i=i 



(this is called dominance ordering in ||5t|; the latter fact follows from Eq. (|T 



and n > n — /X for all A* > 0). This is not at all obvious for the eigenfunctions 
obtained with the novel algorithm (but of course should be true as well, at least 
for non-degenerate eigenfunctions). 

In both algorithms it is important to rule out the occurrence of resonances, but 
the reason for that is different (cf. Eq. (|2T| ) with Lemma |1] above, and observe 
the different sign of /-i). 

In the novel algorithm the restriction in Eq. (|12|) can be dropped, and in fact the 
solutions thus obtained are relevant in the elliptic case ||L1[ . There seems no way 
to drop this restriction in the standard algorithm. 

From Sutherland's algorithm it seems somewhat surprising that the eigenvalues 
all can be written in the simple form En = -Pf > but from the novel algorithm 
this is obvious. 

As discussed in the Introduction, the novel algorithm can be generalized to the 



elliptic case LI, L2 . 
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Appendix A: Proof of Fact 2 



We note that, 



F(x;y) = ^,(^)^,(y) 



(Al) 



with \E'o in Eq- (0) and ip in Eq. (H). Using H'^q = Eq'^q |[Su1|| and the Leibniz rule of 
differentiation we obtain 

i/(x)F(x; y) = + ^ 0(x, - Xk)(^{xj - y,) + 
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j,k j,k,e 



where 



(/)(r) = ip' {r) /ip{r) = |cot(r/2) 
(the prime indicates differentiation). Thus 

"if(x)-if(y)lF(x;y) = A2(-)F(x;y) 



(A2) 



with 

i,£ k^j j,k e^k 

where '[x y]' means 'the same terms but with x and y interchanged' (we used 
that all terms which are even under [x y] cancel). Relabeling indices and using 

(h(r) = —(i>(r) we rewrite 



(■) = ^ ^ ^(p{xj - Xk)(f){xj - ye) + (f){xk - Xj)(f){xk - yi) - 
-(j){xi - yj)(p{xi - ?/fc)] - X ^ y] = ^ ^ ^(f){xj - Xk)^{xj - ye) + 



+(p{xj - Xk)(t>{ye - Xk) + (pive - Xj)(f){ye - Xk) - x^y 
We now can use the trigonometric identity 

cot(x) cot(y) + cot(x) cot(2) + cot{y) cot(2;) = 1 if x + y + z = 0, (A3) 
which shows that 



{xj - Xk)(f){xj - yi) + (j){xj - Xk)(j){ye - Xk) + 0(2/^ - Xj)(j){yt - Xk) 



and thus proves (•) = 0. 



□ 



Appendix B: Other proofs 



B.l Proof of Proposition [I| 

We first observe two simple but useful facts. Firstly, the relation in Eq. (§) remains 
true if we replace F(x; y) by 

F'(x; y) = ce'^^f-i("^-J^^) F(x; y) (Bl) 
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for arbitrary constants P G M and c G C. [To see this, introduce center-of-mass 
coordinates X = Ylj=i^j/^ 

) for i = 2, . . . , N, and similarly for 
y. Then if(x) = —d'^/dX'^ + ifc(x'), and similarly for -ff(y). Invariance of Eq. 
under F e-'^nx-Y)N p ^j^^g follows from {d/dX + d/dY)F{x; y) = 0, and the latter 
is implied by the obvious invariance of -F(x; y) under Xj — > Xj + a, yj yj + a, a E M. 
The invariance of Eq. (j^) under F cF is trivial, of course]. Secondly, the variables 
yj in Eq. (^) need not be real but can be complex. 

As mentioned, we intend to perform a Fourier transformation of the identity in Eq. 
(P), i.e. apply to it (27r)^^/ d^ye'^'^ with suitable momenta P. We need to do this 
with care: firstly, the differential operator H{y) has singularities at points yj = yk, 
and secondly, the function F(x;y) is not periodic in the variables yj but changes by 
phase factors under yj yj + 27r. We therefore need to specify suitable integration 
contours for the y/s avoiding the singular points, and we need to choose the Pj so as 
to compensate the non-periodicity. To do that, we replace the real coordinates yj by 

= Vj - (B2) 

with e > a. regularization parameter: as we will see, we can then integrate along the 
straight lines from yj = —n to vr and after that perform the limit e i 0. Since for all 
j < k, Zj — Zk = yj — yk + ^^kj with Skj = {k — j)e > 0, we can use 

sm[{y + ie)/2] = ie''^/2g-W2+e/2(^ _ ^iy-e^ (33^ 

for e > 0. Taking the log of this identity and differentiating we obtain (l/2cot[(|/ + 
i£:)/2] = — i/[l — exp {iy — e)]. Expanding the r.h.s. in a geometric series and differen- 
tiating once more yields 



4sin^[(^/ + i£)/2] 



J2 i^e^'^y-"'. (B4) 



This accounts for all singularities and branch cuts in a consistent way. To determine 
the suitable P use Eq. (BBf) and compute 



F(x;z) = (---)*o(x)P^(x;y) 



with \E'o(x) in Eqs. 



^^(x; y) = - - ; ^ (B5) 
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a function periodic in all the and 



. N(N~l)/2-N^V\ p-iA(%-yfc)/2+A(fc-j)e/2 



2^ I TjN -i\(x,-yu)/2+\ke/2 

1 lj,fc=l 



2 



= const. e^^^Sf=i(^^-%)/2g-iAEf=i(^+i-2j)i/,/: 

(we used '^^^^iVj — Vk) = + 1 — 2j)yj). We thus see that we can choose P and 

c in Eq. (|B]] ) such that 

F'(x; z) = e-'^^f=i[(^+i)/2-j1y,p.(^. y)^p(x) . (B6) 

We need to choose the Fourier variables P = (Pi, . . . , Pjv) such that e'^'^F'(x; z) is 
periodic in all i/j, and this implies 

p. = + A [i (AT + 1) - j] , G Z. (B7) 

We now can apply (27r)~^ / d^y e^^ ^ to the identity if(x)F'(x; z) = H{z)F'{x; z). 
We recall 

^^^^ ^ " ^ ^ ^ ^ § 4 sin2[(y, - + i(A; - j»/2] 



and use Eq. (|B4|) . After taking the limit £ J, we obtain Eq. (^): the l.h.s. is obvious 
(note that F is the Fourier transform of F'). The r.h.s. has two terms. The first one is 
equal to Ylj PfP comes from the derivative terms which we evaluated by partial 
integration. The second term is obtained from the 1/ sin^-terms in H{z)F' which we 
computed using Eq. (|B4D . □ 



B.2 Proof of Lemma [I| 

We write (n + = rij + fij with fij in Eq. (|TB|). Thus Eqs. ( PHj ) and ( P5| ) imply, 

Er.+^ -E^ = Y, + 2/i,(n, + A[i(Ar + 1) - j])) , 

j 

and with Eq. (H) we obtain Eq. (H). □ 
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B.3 Proof of Lemma g 



It is straightforward to evaluate Vni'x.) in Eq. (|28|) by expanding all terms in Tay- 
lor series (using the binomial series) and then performing the yj integrations which 
corresponds to a projection onto the ?/j-independent terms. The results is 



l<j'<k'<Nj,k=l ^ ^ 

where the sum is over all non-negative integers Hkk' and Uj£ such that 

N j-1 N 

e=i 1=1 i=j+i 



Recalling the definition of in Eq. (|TT]) we obtain Eqs. (plD -fpq). 

We now argue that this latter system of equations can have solutions only if the 
conditions in Eq. (^) all hold, which implies that otherwise Vn is zero. To see this 
we add up the last N + 1 — k relation in Eq. (|B9| ) {k = N, N — 1, . . . ,1), and by a 
relabeling of indices we obtain 

N N N k-l N 

Y'^^ = YY + YY^^'j 

j=k j=k 1=1 i=l j=k 

where the r.h.s. is always manifestly positive. This proves Eq. (pS]). Setting k = 1 and 
comparing with Eq. (|36D we obtain Eq. (^7]). Moreover, for fixed nj, there are at most 
a finite number of different solutions of Eq. (|B9D, implying that Vn is a polynomial. 
To see that we write Eq. ( |B9|) as follows, 

N N j-1 

+ Y ^^'^ = Y + Y (^^°) 

i=j+i e=i 1=1 



and determine possible solutions for decreasing values of j starting at j = N. It is easy 
to prove by induction that there is only a finite number of solutions 

r ^N r ^ ^ j..N^+N(N-l)/2 

i^ijlj,i=i^ wefi<j<e<N ^ ^0 
of this system of equations: For j = N we get 

N N-1 

riN = Y^'^^ ^Y 

1=1 1=1 
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and there is obviously only a finite number of different solutions {i^eN}£=i, {/^w}£=i of 
that equation. If we consider Eq. (|B1(]| ) for some j = Jq < N, the possible solutions for 
{fJ'joAiyjo '^^^^ already determined by the equations for j > jo and, by the induction 
hypothesis, there is only a finite number of them. One thus only has to consider a finite 
number of equations, and each of them obviously has only a finite number of solutions 
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